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Abstract 

We present a method to prove hypergeometric double summation 
identities. Given a hypergeometric term F(n,i,j), we aim to find 
a difference operator L = ao(n)N° + a\{n)N l + • • • + a r (n)N r and 
rational functions R\(n,i, j), R,2(n,i, j) such that LF = Ai(RiF) + 
Aj(i?2^)- Based on simple divisibility considerations, we show that 
the denominators of R± and R2 must possess certain factors which 
can be computed from F(n,i,j). Using these factors as estimates, we 
may find the numerators of R± and R2 by guessing the upper bounds 
of the degrees and solving systems of linear equations. Our method is 
valid for the Andrews-Paule identity, Carlitz's identities, the Apery- 
Schmidt-Strehl identity, the Graham-Knuth-Patashnik identity, and 
the Petkovsek-Wilf-Zeilberger identity. 
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1. Introduction 

This paper is concerned with double summations of hypergeometric terms 
F(n, A function F(n, fa, . . . , k m ) is called a hypergeometric term if the 
quotients 

F(n + 1, k u . . . , km) F(n,ki + l,...,k m ) F(n, k x , . . . , k m + 1) 



F(n, fa, . . . , k m ) F(n,ki,...,km) F(n, fa, . . . , k m ) 

are rational functions of n,fa, . . . , k m . Throughout the paper, we use N to 
denote the shift operator with respect to the variable n, given by NF(n) = 
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F(n + 1) and use A x to denote the difference operator with respect to the 
variable x, given by A X F = F(x + 1) — F(x). For polynomials a and b, we 
denote by gcd(a, b) their monic greatest common divisor. When we express 
a rational function as a quotient p/q, we always assume that p and q are 
relatively prime unless it is explicitly stated otherwise. 

Zeilberger's algorithm [14,17,22], also known as the method of creative 
telescoping, is devised for proving hypergeometric identities of the form 

£>(n, *) = /("), (1-1) 
k 

where F(n,k) is a hypergeometric term and f(n) is a given function. This 
algorithm has been used to deal with multiple sums by Wilf and Zeilberger 
[21]. Given a hypergeometric term F(n, k±, . . . , k m ), the approach of Wilf 
and Zeilberger tries to find a linear difference operator L with coefficients 
being polynomials in n 

L = a (n)N° + o^nJJV 1 + • • • + a r (n)N r 

and rational functions Ri, . . . , R m of n, ki, . . . , k m such that 

m 

LF = '£,A kl (R i F). (1.2) 
i=i 

As noted by K. Wegschaider [20], when the boundary conditions are ad- 
missible, Equation (jl.2j) leads to a homogenous recursion for the multi- 
summations: 

L ^ F(n,k 1 ,...,k m ) = 0. 

When m — 1, L and R\ can be solved by Gosper's algorithm [13, 17]. 
S.A. Abramov, K.O. Geddes and H.Q. Le also provided a lower bound for the 
order r [2,3] and found a faster algorithm [4] compared with Zeilberger's algo- 
rithm. For a survey on recent developments, see [1]. For m > 2, constructing 
the denominators of R±, . . . , R m for the Wilf- Zeilberger approach remains an 
open problem. In a recent paper [16], M. Mohammed and D. Zeilberger used 
the denominator of LF/F as estimates of the denominators of Ri. In an al- 
ternative approach, Wegschaider generalized Sister Celine's technique [20] to 
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multiple summations, and proved many double summation identities. A dif- 
ferent approach has been proposed by F. Chyzak [11,12] by finding recursions 
of the summation iteratively starting from the inner sum. C. Schneider [18] 
presented the Chyzak method from the point of view of Karr's difference field 
theory. 

To give a sketch of our approach, we first consider Gosper's algorithm 
for bivariate hypergeometric terms. Suppose that F(i,j) is a hypergeometric 
term and Pi/qx,p 2 /q 2 are rational functions such that 

F(i,j) = A, ( eMjPft,!) + a, ( . 

We show that under certain hypotheses (Section 2, (H1)-(H3)), the denom- 
inators q\, q2 can be written in the form 

qi(i,j) = v x (i) v 2 (j) v 3 (i + j)v A (i,j)ui(j)u 2 (i,j), 

(1.3) 

= vi(i)v 2 (j) v z (i + j)v 4 (i, j)wx(i) w 2 (i, j), 

such that vi,v 2 , t>4 and u 2 , w 2 are bounded in the sense that they are factors 
of certain polynomials which can be computed for a given F(i,j), see The- 
orem 12.11 Then we apply these estimates to the telescoping algorithm for 
double summations. Suppose that 

LF(n,i,j) = A i (R 1 (n,i,j)F(n,i,j)) + A j (R 2 (n,i,j)F(n,i,j)), 

where 

p /_ _ 1 hi^hj) p / • « 1 f2(n,i,j) 



d(n,i,j) gi(n,i,j) " ' d(n,i,j) g 2 (n,i,j) 

and d(n,i,j) is the denominator of LF(n,i,j)/F(n,i,j). We may deduce 
that gi, g 2 can be factored in the form of (jl.3|) such that Vi,v 2 , v 4 and u 2 , w 2 
are bounded, see Theorem 13.11 Although we do not have the universal de- 
nominators, these bounds can be used to give estimates of the denominators 
gi and g 2 . Then by further guessing the bounds of the degrees of the numer- 
ators of Ri and R 2 , we get the desired difference operator if we are lucky. 

Indeed, our approach works quite efficiently for many identities such as 
the Andrews-Paule identity, Carlitz's identities, the Apery-Schmidt-Strehl 
identity, the Graham-Knuth-Patashnik identity, and the Petkovsek-Wilf-Zeilberger 
identity. 
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2. Denominators in Bivariate Gosper's Algorithm 



For a given bivariate hypergeometric term F(i,j), we give estimates of the 
denominators of the rational functions Ri(i, j), R 2 (i, j) satisfying 



F(iJ) = A i (R 1 (i,j)F(i,j))+A j (Il a (i,j)F(i,j)). (2.1) 



Let 



RlihJ) = / • ■ \ ; R 2 (l,f 



F{i + l,j) r x (i,j) F{i,j + 1) r 2 (i,j) 



(2.2) 



F(i,j) si(i,j)' F(i,j) s 2 (i,j)' 

Dividing F(i,j) on both sides of (j2.1|) and substituting (|2.2j) into it, we derive 
that 

1 _ h(i + j) | r 2 (i, j) f 2 (i, j + 1) / 2 (»,j) , 23 s 

si(« ) j)fl , i(* + 1 )i) ^(^i) s 2 (i, j) & (»\ J + 1) 92{i,j)' 

Let 

=g c d(si(i,j),s 2 (i,j)), f(i,j) =gcd(g 1 (i,j),g 2 (i,j)), 



and 

s'i(^i) = si(i,j)/u(i,j), s' 2 (ij) = s 2 (i,j)/u(i,j) 
9i(i,j) = 9i(iJ)/v(i,j), g' 2 (i,j) = 92(i,j)/v(i,j) 



(2.4) 



We find that in many cases we can restrict our attention to those Ri, R 2 
whose denominators gi,g 2 satisfy the following three hypotheses. We see 
that in the proof of the following theorem, these hypotheses enable us to 
cancel out unknown factors from the multiples of g\ and g 2 so that we can 
obtain an upper bound of g\ and g 2 . Thus, these hypotheses come naturally 
from the requirement of simple divisibility properties. Moreover, it turns 
out that these divisibility requirements are sufficient in many cases to give 
good estimates for the denominators g\ and g 2 . The three hypotheses are as 
follows: 

(HI) Suppose p(i,j) and p(i+ hi, j + h 2 ) are both irreducible factors of g\ (i,j) 
(g 2 (i,j), respectively) for some hi, h 2 G { — 1, 0, 1}. Then they must be 
coincide. 



4 



(H2) gcd(g' 1 (iJ),v(hj))=Scd(g' 2 (iJ)Mhj)) = 1 - 
(H3) For any integers hi, h 2 G {—1, 0, 1}, 

gcd(g[(i + h u j + h 2 ),g' 2 (i,j)) = 1. 

For example, the following functions satisfy the above hypotheses: 

9l (i,j) = (2n-2i+l)(n-i+l)(j + l) 2 , g 2 (i,j) = (2n-2i+l){n-i+l){i+l) 2 . 

Remarks. 

1. Hypothesis (HI) looks like requiring that g\ and g 2 are shift-free (see 
Abramov and Petkovsek [5]). However, only the shifts of ±1 are consid- 
ered and shift invariant factors are admissible. For example, we allow 
that gi(i,j) = (i + l)(i + 3) or g x {i,j) = i+j. 

2. According to [6], gcd(^i(i, j), g x (i + h u j + h 2 )) and gcd(g 2 (i, j), g 2 (i + 
hi,j + h 2 )) can factor into integer-linear factors for hi,h 2 being not 
both zero. 

3. Hypothesis (H3) is to require that g[/ g 2 are shift-reduced (see also [5]) 
respect to the shifts of ±1. 

Under the above hypotheses, we have 

Theorem 2.1 The denominators g±(i, j), g 2 (i, j) can be factored into poly- 
nomials: 

gi(ij) = V!(i)v 2 (j)v 3 (i + j)v 4 (i, j)m(j)u 2 (i, j), 
92(i,j) = V!(i)v 2 (j)v 3 (i + j)v 4 (i, j)wi(i)w 2 (i,j), 

such that 



vi(i) \ r 1 (i - l,j)s 2 (i - 1, j), 




(2.5) 


v 2 (j) \r 2 (i,j - l)s[(i,j - 1), 




(2.6) 


v 4 (i,j) | gcd (n(i - l,j)s' 2 (i - l,j),r 2 (i,j - 


l)s[(i,j-l)), 


(2.7) 


u 2 (i,j) | gcd (s l (i,j)s 2 (i,j),r l (i - l,j)s 2 (i - 


-u)), 


(2.8) 


w 2 (i,j) I gcd (s 2 (i,j)s[(i, j),r 2 (i, j - l)s[(i,j - 1)). 


(2.9) 
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Proof. Substituting (JUJ) into O, we get 

, ri(i,j) + 



s'i(i,j)u(i,j) g[(i + 1, j)u(z + 1, j) g[(i, j)v(i,j) 
r 2 (ij) /2(i,j + l) /2(*,j) 



s' 2 (h j) g' 2 (i, j + l)v(i,j + 1) g' 2 (i,j)v(i,j)' 

That is, 

si(ij)s' 2 (ij)gi(hj)g' 2 (hj)gi(i + l,j)g 2 (i,j + 1) 
= + l,j)ri(i,j)s 2 (i,j)gi(i,j)g 2 (i,j)g 2 (i,j + 1) 
-fi(hj)si(i,j)s' 2 (i,j)g' 2 (i,j)g 1 (i + l,j)g 2 (i,j + 1) 

+/2(», j + l)r 2 (i, j)gi(i,j)g' 2 (hj)gi(i + 1 J) 

-f2(i,j)si(i,j)s' 2 (i,j)g' 1 (ij)g l (i + l,j)g 2 (i,j + 1). 

1. Suppose that p(i,j) is an irreducible factor of v(i,j), and for some non- 
negative integer Z, 1 1>. Note that p(i + h±,j + h 2 ) is also irreducible. 
Since 

gcd(p(i + l,j),A(* + l,j')) = gcd(p(i,j + l),/2(i,J + l)) = 1, 
we have 

p'(z + 1, j) \ri(i,j)s 2 (i,j)g 1 (i,j)g' 2 (i,j)g 2 (i,j + 1) 

and 

j + 1) | r 2 (i,j)s' l (i,j)g l (i,j)g' 2 (i,j)g l (i + 

There are three cases: 

• p(i,j) is a polynomial depending only on z. Then gcd(p(z + 
1) j)i Qi(h j)) — 1- Otherwise, by hypothesis (HI) we have that 
p(i + 1, j) = p(i,j) is independent of z, which is a contradiction. 
Similarly, gcd(p(i+l, j), g 2 (i,j)) = 1. Since p(i,j) is a polynomial 
depending only on z, we have 

gcd(p(i + l,j),(/ 2 (i,j + l)) = gcd(p(z + 1, j + l),g 2 (i, j + 1)) = I. 
Hence, 

p l (i + 1, j) |n(z, j)s' 2 (i,j). 

Let fi(z') denote the product of all irreducible factors of v(i,j) 
that depend only on z. Then we have (J2.5|) . 
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• p(i,j) is a polynomial depending only on j. The same discussion 
leads to 

p l (i,j + 1) | r 2 (i,j)s[(ij). 

Let v 2 (j) denote the product of all irreducible factors of v(i,j) 
that depend only on j. Then we have (J2.6J) . 

• p(i,j) is a polynomial depending both on i and on j. Then either 

p(i + l,j)=p(i,j + l) (2.10) 

or 

gcd(p(i+l,j),p(i,j + l)) = l. (2.11) 

In the former case, p(i,j) is a polynomial of % +j (see [6, Lemma 
3] or [15, Lemma 3.3]). For this case we do not have a bound. We 
denote by v 3 (i + j) the product of all irreducible factors p(i,j) of 
v(i,j) that satisfy (|2.1(J|) . In the later case, by hypothesis (HI), 
we have 

gcd(p(z + l,j),g 1 (i,j)g' 2 (i,j)g 2 (i,j + 1)) = 1 

and 

gcd(p(i,j + l),g 1 (i,j)g' 2 (i,j)gi(i+lj)) = 1. 

Thus, 

P l {hj) I gcd (n(i - l,j)s' 2 (i - l,j),r 2 (i,j - l)s[(i,j - 1)). 

Let V4(i,j) denote the product of all irreducible factors p(i,j) of 
v(i,j) that satisfy ([2. lip . Then we have (|2.7p . 

2. Suppose p is an irreducible factor of g[ and p \g[ for some non-negative 
integer I. If p(i,j) \ v(i,j + 1), then p(i,j — 1) | v(i,j). By hypothesis 
(HI), p(i,j — 1) = p(i,j), which implies p(i,j) \ v(i,j), contradicting 
to hypothesis (H2). Noting further that by hypothesis (H3), for any 
h u h 2 G {-1,0,1}, 

gcd(f!(i,j),gi(i,j)) = gcd(g[(i,j),g' 2 (i + h u j + h 2 )) = 1, 

we have that 

P l (hj) I s 1 (i,j)s' 2 (i,j)g 1 (i + l,j). 
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If p(i + l,j) \v(i,j + 1), then by hypothesis (HI), p(i + l,j — 1) = 
p(i,j), which implies p(i,j) \ v(i,j), contradicting to hypothesis (H2). 
Therefore, by hypothesis (H3), 

p\i + 1, j) | r 1 {i,j)s' 2 (i,j)g 1 (i,j). 

There are two cases: 

• Pihj) — P(i + 1) j)- Then p(i,j) is a polynomial depending only 
on j. For this case we also do not have a bound. We denote by 
Ui(j) the product of all irreducible factors of g[(i,j) that depend 
only on j. 

• gcd(p(i, j),p(i + 1, j)) = 1- Then by hypothesis (HI), 

gcd(p(i,j),g x (i + l,j)) = gcd(p(i + l,j),g x (i,j)) = 1, 
and hence, 

P l (i,j) I gcd (si(z, j)s' 2 (i,j),n(i - l,j)s' 2 (i - 1, j)). 

Let U2(i,j) denote the product of all irreducible factors p(i,j) of 
g[(i,j) such that gcd(p(i, j),p(i + 1, j)) = 1. Then we have ()2.8|) . 

3. Similarly, suppose p is an irreducible factor of g' 2 and p l \g' 2 for some 
non-negative integer /. Then either p(i,j) is a polynomial depending 
only on i or 

P l (hj) I gcd (s 2 (i,j)s[(i,j),r 2 (i,j - l)s[(i,j - 1)). 

Let Wi(i) denote product of irreducible factors of g' 2 (i,j) that depend 
only on i and w 2 (i,j) denote the product of the rest irreducible factors 
of g' 2 (i,j)- Then we have (j2.9J) . I 

Note that u 2 (i,j) have no factors which are free of i and w 2 (i,j) have no 
factors which are free of j. We will need this property later for the algorithm 
EstDen. 
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3. Denominators in Our Telescoping Method 

We are now ready to estimate the denominators of R\ and R 2 in our tele- 
scoping method. 

As in the case of single summations, the telescoping algorithm for double 
summations tries to find an operator 

L = cio(n) + a\{n)N + • • • a r (n)N r 

and rational functions Ri(n,i, j), R 2 (n,i, j) such that 

LF(n, i,j) = Ai(i2! (n, i, j)F(n, 2, j)) + A j (R 2 (n, 2, j)F(n } i, j)). (3.1) 



(3.2) 



Let 

F(n,i + 1, j) _ rx(n,z, j) F(n,i,j + 1) _ r 2 (n,i,j) 
F(n,i,j) si(n,i,jY F(n,i,j) s 2 (n,i,j)' 

and d(n,i,j) be the common denominator of 

F(n + F(n + r,i,j) 

F(n,i,j) ' '"' F(n,i,j) 

Then there exists a polynomial c(n,i,j), not necessarily being coprime to d, 
such that 

LF(n,i,j) F(n + l,i,j) _ c(n,i,j) 

F(n,i,j) F ^ 1 ^) dfaiJY 

Note that c is related to the polynomials CLq, Gil, . . . , O r but g? is independent 
of them. 

Now, (|3.1j) can be written in the form of (|2.1jl : 
LF(n, i,j) = Ai(R[ (n, 2, j)LF(n, 2, j)) + A,(i^(n, 2, j)LF(n, 2, j)), 
where 

Ri{n,i,3) = Ri{n,i,j)— — — - and R 2 {n,i, j) = R 2 (n,i, j] 



c[n,i,j) ' cln,t,j) 
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This suggests us to assume 

Ri\n,i,j) = T — — — r and R 2 (n,t,j) = — — — - — — 

d{n,i,j) gi{n,i,j) d{n,i,j) g 2 (n,i,j) 

(3.4) 

where f\ and g\ (f 2 and g 2 , respectively) are relatively prime polynomials. 

Since the following discussion is independent of n, we omit the variable n 
for convenience. For example, we write Ri(i,j) instead of Ri(n,i,j). Using 
these notations, we have 

Theorem 3.1 Suppose the polynomials g\ and g 2 in (|3.4j) satisfy the hy- 
potheses (H1)-(H3) . Suppose further that for any hi, h 2 £ { — 1, 0, 1}, 

gcd(gi(i,j),d(i + h u j + h 2 )) = gcd(g 2 (i,j),d(i + h u j + h 2 )) = 1. (3.5) 

Then gi(i, j), g 2 (i, j) can be factored into polynomials: 

9i(i,j) = vi{i)v 2 (j)v^(i + j)vSd)ui(])u 2 (i, j), 
g 2 (hj) = Vx(i)v 2 (J)v s (i + j)v A (i, j)w 1 (i)w 2 (i, j), 

such that 

v^i) | n(i - l,j)s' 2 (i - l,j), 
v 2 (j)\r 2 (i,j - l)s[(i,j - 1), 

Mhj) I gcd (n(i - lj)s' 2 (i - l,j),r 2 (i,j - l)s[(i,j - 1)), 
MhJ) I § cd { s i(hj) s ' 2 (hj), r i(i - l,j)s' 2 (i - 1, j)), 
w 2 (i,j) | gcd (s 2 (i,j)s[(i,j),r 2 (i,j - l)s[(i,j - 1)), 

where 



s'i(hj) = si(i,j)/gcd(s 1 (ij) 1 s 2 (i,j)), 
s'Sd) = s 2 (i, j)/gcd(si(i, j),s 2 (i, j))- 



(3.6) 



Proof Substituting (|3.4|) into (|3.1|) and dividing F(i,j) on both sides, we 
obtain 



c{hi) _ /i(z + i,j) fi(i,j) _ 

j) 

(3.7) 



d{i,j) s 1 (i,j)d(i + l,j)g 1 (i + l,j) d{i,j)g 1 (i,j) 
, r 2 (i,j) f 2 (i,j + 1) h(i,j) 



s 2 (i,j) d(i,j + l)g 2 (i,j + 1) d(i,j)g 2 (i,jy 
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i.e., 



c / i .x = r 1 (i,j)d( y i,j) f\{i + l,j) _ fi(i,j) 
si{i,j)d(i + l,j) gi(i + l,j) gi(i,j) 
r2(i,j)d(i,j) f 2 (i,j + l) _ f 2 (i,j) 
s 2 (i,j)d(i,j + l)g 2 (i,j + 1) #2(2, j)' 

Let 

^i(^i) = n(i,j)d(i,j), hihj) = s 1 (i,j)d(i + 1, j), 
r 2 (i,j) = r 2 (i,j)d(i,j), s 2 (i,j) = s 2 (i,j)d(i,j + 1). 

All discussions in the proof of Theorem 12.11 still hold. Thus, we have 

Vx{i) |fi(i - l,j)s' 2 (i - 1, j), 
v 2 (j) I f 2 (i,j - j - 1), 

^4(2, j) I gcd (f x (z - l,j)s' 2 (i - l,j),r 2 (i,j - l)s[(i,j - 1)), (3.8) 
u 2 (i,j) I gcd (s 1 (i,j)s' 2 (i,j),f 1 (i - l,j)s' 2 (i - 1, j)), 
w 2 (i,j) I gcd (5 2 («,j) 5 'i(^j)^2(^J - j - 1)), 

where 

&i(iJ) = si(i,j)/gcd(s 1 (i,j),s 2 (hj)), 
s' 2 (i,j) = s 2 (i, j)/gcd(§i(i, j),s 2 (i, j)). 

Since we have (|3.5j) . we may replace ri,§i,f 2 ,S2 by ri,Si,r 2 ,S2 in ()3.8j) . re- 
spectively. I 



4. A Telescoping Method for Bivariate Hypergeometric Terms 

Theorem 13 . II enables us to choose the denominators in the telescoping algo- 
rithm. Basically, we will use certain factors appearing in the bounds of the 
denominators as estimates of the denominators. In many cases, this approach 
seems to work quite efficiently although we are not able to give a formula 
to bound the denominators because certain factors are not bounded in The- 
orem 13.11 Roughly speaking, the divisibility considerations in our method 
serve as a guide to guess the factors in the denominators. In fact, the esti- 
mated denominators are much smaller than the theoretical bounds given by 
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Theorem 13.11 Only u 2 (i,j) and w 2 (i,j) are set to their theoretical bounds, 
while v 2 (j),v 3 (i + j),t> 4 (i, j) are set to 1, U\{j) and Wi(i) are set to factors 
of si(i, j)s' 2 (i, j), and vi{i) is set to a factor of its theoretical bound. See the 
following algorithm EstDen. 

Algorithm EstDen 

Input: A hypergeometric term F(n,i,j). 

Output: Estimated denominators gi(i,j) and g 2 (i,j) for bivariate Gosper's 
algorithm. 

1. Calculate r±, r 2 , Si, s 2 , s i> 

s' 2 defined by (jSH and (EHD) : 

2. Set 

^i(i) := the maximal factor of ri(i, j)s' 2 (i, j) depending only on i\ 
v 2{j) '■= the maximal factor of r 2 (i, j)s[(i, j) depending only on j; 
and 

v(i) := gcd(wi(i - l),v 2 (i - 1)); 

3. Set 

U\(j) := the maximal factor of Si(i, j)s 2 (i, j) depending only on j; 
Wi(i) := the maximal factor of Si(i, j)s 2 (i, j) depending only on i; 

4. Set u 2 (i,j) to be the maximal factor of 

gcd(si(i, j)s' 2 (i, j),n(i - l,j)s' 2 (i - 1, j)) 

which depends on i; 

Set w 2 (i,j) to be the maximal factor of 

gcd(s 1 (i,j)s^ 2 (i, j),r 2 (i, j - l)s[(i,j - 1)) 

which depends on j. 

5. Return gi(i,j) := v(i)ui(j)u 2 (i, j) and g 2 (i,j) := v (i)wi(i)w 2 (i, j). 

Remark. Let f(i,j) be a polynomial in i, j and a be a new variable. Then 
the maximal factor of f(i,j) depending only on i can be obtained by 

gcd(f{i,j),f(i,j + a)), 
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and the maximal factor of f{i,j) depending on i can be obtained by 

f(i,j)/ gcd(f(i,j),f(i + a,j)). 



We are now ready to describe our telescoping method for double summa- 
tions: 

Method BiZeil 

Input: A hypergeometric term F(n,i,j). 

Output: An operator L and rational functions Ri and R 2 such that (|3.1j) 
holds if the algorithm succeeds. 

1. Using algorithm EstDen to obtain g\ and g 2 - 

2. Set the order r of the linear difference operator L to be zero. 

3. For the order r, calculate the common denominator d(n,i,j) of 



(If r = 0, then take d(n,i,j) = 1.) 

4. Set the degrees of f\ and f 2 to be one more than those of d • g\ and 
d ■ g2, respectively. 

5. Solve the equation ()3.7|) by the method of undeterminate coefficients 
to obtain ao, ai, . . . , a r and /1, f 2 . 

6. If ai 7^ for some % G {0, . . . , r}, then return L, fi/(d ■ g±), fij [d ■ g%) 
and we are done. 

If dj = for all i 6 {0, . . . , r}, but deg f\ — deg(d- gi) < 2, then increase 
the degrees of fi and fi by one and repeat step 5. 

Otherwise, set r := r + 1 and repeat the process from step 3. 



1. In many cases, gi{i,j) and g^ihj) can be further reduced by cancelling 
a factor of degree 1 and a factor of degree 2 from g x and g 2 , respec- 
tively. In our implementation we first choose two arbitrary factors and 



F(n + l,i,j) 
F(n,i,j) 



F(n + r,i,j) 



F(n,i,j) 



Remarks. 
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use the reduced g\ and g 2 . When it fails, we then try the unreduced 
ones. This cancellation may reduce the time of calculation if we are 
lucky For example, for the Andrews-Paule identity (see Example 1), 
the estimated denominators given by Theorem 13.11 by algorithm Est- 
Den, and by reduction are, respectively, 



gi (i,j) = (2n - 2i + l)(n - i+ l)(2n - 2j + l)(n-j + 1)(* + j) 2 (j + l) 2 , 
g 2 (i,j) = (2n - 2t + l)(n - i + l)(2n - 2j + l)(n - j + l)(t +i) 2 (i + l) 2 ; 

9l (i,j) = (2n - 2i + l)(n - i + l)(j + l) 2 , 
52(t,i) = (2n - 2t + l)(n - t + l)(i + l) 2 ; 



^(i, j) = (2n - 2z + l)(j + l) 2 , g 2 (i,j) = (2n - 2i + l)(n + 



The calculation times are 116 seconds, 5 seconds and 0.6 second, re- 
spectively. We should note that since our method is heuristic and it 
applies only to particular cases, we are more interested in the compu- 
tation results which are verifiable. So we cannot claim the efficiency of 
the method or its applicability. 

2. In all the following examples except Example 4, the degree of the nu- 
merator of Ri is one more than that of the denominator. While 
in Example 4, the difference is two. 

The degree bounds can be interpreted as follows. Let t\, t 2 , t$, £4 be the 
four terms of the right hand side of (|3.7J) after multiplying the common 
denominator. In most cases, the leading terms of £1 and t 2 (£3 and £4, 
respectively) are cancelled. 

3. There is a way to speed up the computation in Step 5. Given g 1 and 
g 2 , we may derive part of the factors of fi and f 2 by divisibility. For 
example, suppose (|3.7|) becomes 



after substituting and simplification. Suppose further that D(i,j) is 
the common denominator of the above equation. Then we immediately 
have that fi-D/wi is divisible by q x = gcd(cD/d, u±D/vi, u 2 D/v 2 , D/w 2 ) 



and 



c(i,j) = ui(i,j) 
d(i,j) v^ij) 




w 2 {i,j 



f2(i,j) 
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and fi(i+l,j)-UiD/vi is divisible by q 2 = gcd(cD/d, D/wi, u 2 D/v 2 , D/w 2 ), 
and hence, 

qi , <?2 

and 



gcd(D/wi,gi) gcd{uiD/v 1 ,q 2 ) 
are factors of fi(i,j) and /i(z + 1, j), respectively. 



5. Examples 

In the following examples, let F denote the summand of the left hand side 
of the identity. 

Example 1. The Andrews- Paule identity: 

i=0 j=0 



It was proved by G. Andrews and P. Paule [7, 8] by establishing a more 
general identity 

LfJ LfJ x . , .\ 2 



EE 



'i+ j\ 2 fm + n-2i-2j\ [^^\\[^^\\ 

l=0 , =0 v*A n-2i J ~ TflWlWIW' 

Using the method BiZeil, we can deal with ()5.1j) directly. In fact, we have 

9l {iJ) = (2n-2m)(n-*+l)(j + l) 2 , g 2 {i,j) = {2n-2i+l){n-i+l){i+lf . 

Cancelling the factors (n — i + 1) and (i + l) 2 from gi(i,j) and g 2 (i,j), 
respectively, we obtain 

gxiij) = (2n-2« + l)(j + l) 2 and g 2 (i,j) = (2n - 2z + l)(n - z + 1). 

Finally, we get 

(2n + z, j) = A^F^n, + AjR 2 F{n, i, j), 

where 

i 2 (6n 2 + 5n + 1 + 6jn 2 + — j — in + 2m 2 — 2z — 4j 2 n — 2j 2 — 3ij — 4zjn) 



(2n -2i + 1)(1 + j) 2 
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— 2n 2 + 2jn 2 + 6m 2 + 9in + 3jn — Aijn — Ai 2 n — n + j — 3ij + 2i — Ai 2 
2 ~ (2n - 2i + 1) ' 

which are the same as given in [20, p. 85]. Summing z, j = 0, . . . , n, we get 

n n 

(272 + l)^^F(n,*,j) 

n n 

= {R2F{n, i,n + l)- R 2 F{n, i, 0)) + (RxF(n, n + l,j)- R\F{n, 0, j)) 

i=0 j=0 

Note that there is only one nonzero term RiF(n,n + l,n) of the second 
summation. While applying Gosper's algorithm to the first summand, we 
obtain 

n 

^2 (R 2 F{n, i,n + l)- R 2 F(n, i, 0)) = G(n + 1) - G(0). 

i=0 

where 

(-2n + i - 1) (-An + 2i- l)i (An - 2i N 
^ ~ -1 + 2% - 2n \2n - 2z, 

Simplifying G(n + 1) - G(0) + RiF(n, n + 1, n), we finally get (fO) . 
Example 2. Carlitz's identity [10] (see Also [21, Example 6.1.2]): 



i=0 j=0 v J/ z=o 



We have 

= (j + l) 2 (-n + j), <7 2 (i,j) = (« + l) 2 (-n + 0. 

Cancelling the factors (— n + j) and (i + 1)(— n + z), we obtain 

<7i(«,j) = (j + l) 2 and g 2 (i,j) =i + l. 

Notice that the common denominator of 

F(n + l,i,j) F(n + 2,i,j) 

TTi \ — an d ^7 s — 

F{n,i,j) F{n,i,j) 
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is (— n + i — 1 + j) 2 (—n + i — 2 + j) 2 . We finally get 

L = (in + 6) - (8 + 5n)N + (n + 2)N 2 , 

and 

R x = ( - i 2 (- n + z — 1)(36 — 10ji 2 n - 13j 2 m + 60j 2 + 60ji - 2i 2 - 38j 2 i 

- 8ji 2 + 10i 3 + 3Qn 3 - 11m 3 - Ujn 3 - 2z 4 - 92jn 2 + 8i 2 n - 80m + 5j V 
+ 8j 2 i 2 + 88jm + 42j 2 n - 172 jn + 24jin 2 + 5i 2 n 2 + 3i 3 n - 54m 2 + 88n 2 + Aj 3 n 

-90j + 6j 3 -40i + 5n 4 + 90n)) / ((-n + i - 1 + j) 2 (-n + i -2 + j) 2 (j + 1) 2 ) , 

R 2 = ((64 - 19 ji 2 n - 6j 2 ni + Uj 2 + 74ji + 54i 2 - 10j 2 i - 36ji 2 + 2i 3 + 39n 3 

- 16m 3 - 9jn 3 - Ai 4 + Qji 3 - 53 jn 2 + 50i 2 n - 176m + 4j 2 n 2 + Aj 2 i 2 + 5n 4 
+ 83jm + 16j 2 n - lOOjn + 22 jm 2 + llz 2 ra 2 + Ai 3 n - 93m 2 + 112n 2 - 60 j 

- 108z + 140n)(-n - 1 + j)) j ((-n + % - 2 + jf{-n + i - 1 + j) 2 ), 



such that 

LF(n,i,j) = Aji?iF(n, i, j) + AjR 2 F(n, i, j), (5.2) 

By summing ()5.2j) over i,j from to n, one derives that L annihilates 
the double sum on the left hand side. It is easily seen that the right hand 
side can be annihilated by f(n + 2)iV — (4n + 6))(iV — 1), which is exactly 
L. Then the identity follows from the initial values n — 0, 1. 

The proofs of the following examples are similar to that of Example 2. We 
only need to give #1,(72 and L,Rx,R 2 . Then these identities can be verified 
by checking the initial values. 

Example 3. Carlitz's identity [9] (see also [21, Example 6.1.3]): 



EE 

i=0 j=0 



i + j\ / m — i + j\ fn — j + i\ (m + n — i—j 
i )\ j )\ i )\ m-i 



(m + n + 1)! 1 fm\ fn 



m\n\ ~^ 2k + 1 V k J \k 

k 



By cancelling the factors (1 + j) and (i + l) 2 , we obtain 

9i(hj) = (n-j+i)(l+j) and g 2 (i, j) = m - i + j 
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Notice that the common denominator of 

F(n + l,i,j) F(n + 2,i,j) 

F(n,i,j) " F(n,i,j) 

is (—n + j — 1) 2 (— n + j — 2) 2 , which is denoted by d(i,j). We finally get 

L = 2{m + 3 + n)(2 + m + nf- 

(3m + 2nm + An 2 + 14 + 15n)(n + m + 3)iV + (2n + 5)(n + 2) 2 N 2 , 

and the denominators of R\,R2 are d(i, j)g\(i, j) and d(i,j)g2(i,j), respec- 
tively. The degrees of denominators and numerators of R±,R2 are both less 
than those given in [21]. 

Example 4. The Apery-Schmidt-Strehl identity [19]: 

By cancelling the factors (— j — 1 + and (i + l) 2 , we obtain 

gi(h j) = (~j ~ 1 + i) 2 and g 2 (i,j)=i + l 
Notice that the common denominator of 

F(n + l,i,j) , F(n + 2,i,j) 

— — ; r ctllQ — — ; r 

F(n,t,j) F(n,i,j) 
is (n + 2 — j)(n + 1 — j). We finally get 

L = (n + l) 3 - (3 + 2n)(17n 2 + 51n + 39)A^ + (n + 2) 3 iV 2 , 

and 

i?i = (-2i 2 (3 + 2n)(-10 + 30j 2 -49n 2 -j 3 -4n 4 -24n 3 -2n¥ + n 2 i-6m 2 
+ 3m + 3nji + n 2 ji + 3j 2 i 2 - 3j 3 i + 3ji - Ai 2 - 2j 2 i - 2ji 2 + lln 2 j 2 + 6n 2 j 

+ 33nj 2 + 18nj-Qf + 2t + 15j-39n)) I (( n + 2-j)(n+l-j)(-j -1+i) 2 ), 
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R 2 = (2(-j + z)(3 + 2n)(-8n 2 i - An 2 i 2 - An 2 ji + 4n 2 j + 4n 2 j 2 + 12nj 
- 12nji - 24m + Ylnf - Ylni 2 + 12j 2 - 4jz 2 + f + Qj 2 i 2 - 3j 4 + 8j 

+ 5j 2 i - 8i 2 + 3fi - m - 16ji)) /((n + 2 - j)(n + 1 - + 1)). 

The rational functions R±,R2 are simpler than those given in [19]. The 
operator L was used by Apery in his proof of the irrationality of ((3) and 
Chyzak and Salvy obtained it using Ore algebras [12]. 

Example 5. The Strehl identity [19]: 

By cancelling the factor (— 3i — 3 + 3z — 2 + j) from #2, we obtain 

9i(i,j) = (j + l-*) 3 and ^ 2 (?,j) = (-3«-l+j)(^ + l) 3 - 
Notice that the common denominator of 

F(n + l,i,j) F(n + 6,i,j) 

F(n,i,j) ' '"' F(n,i,j) 

is (n + 1— j)(n + 2— j) - • • (n + 6 — j), which is denoted by d(i,j). We finally 
get a linear difference operator L of order 6 and the denominators of R\ , R2 
are d(i, j)gi(i, j) and d(i, j)g2(i, j), respectively. The operator L is the same 
as the operator obtained by applying Zeilberger's algorithm to the right hand 
side of Q. 

Example 6. The Graham-Knuth-Patashnik identity [14, p. 172]: 

EE(-)-(i + + :) c + r_r *) - <-o U - - - 

j A; 

(5.4) 

By cancelling the factor (j + l)(j + 1 — I) from 5f 2; we obtain 
g 1 (j,k) = (k + l)(k + l + l) and g 2 (j,k) = l. 

Notice that the denominator of ^~fe~yf) is r — j + 1 , which is denoted by 
/c). We finally get a linear difference operator with respect to the variable 

r: 

L = (r + n + l)(n + s + / — m — r) + (r — / + l)(r — s)i? 
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and the denominators of Ri, R 2 are d(j, k)g\(j, k) and d(j, k)g 2 (j, k), respec- 
tively. Then (|5.4j) follows from the evaluation of the initial value (r = 0) by 
Zeilberger's algorithm: 

(fc + z)(fc)( m ) = ^ (n + l)(m-n-l)' 
Example 7. The Petkovsek-Wilf-Zeilberger identity [17, p. 33]: 




(5.5) 

By cancelling the factors s + 1 and (r + l) 2 , we obtain 

gi(r, s) = (n + r)(n + 1 — r)(s + 1) and g2,(r, s) = (n + r)(n + 1 — r). 

Notice that the common denominator of 

F(n+l,r,8) , F(n + 2,r,s) 

77? \ — an d 777 \ — 

F{n,r,s) F{n,r,s) 

is 

(n+l)(n+2)(n+l-r)(n+2-r)(n+l-s)(n+2-s)(n-r-s+l)(n+2-r-s), 

which is denoted by d(r, s). We finally get 

L = 4(4:n + 5)(4n + 3)(n + 1) + 2(2n + 3)(3n 2 + 9n + 7)N - (n + 2) 3 N 2 

is a linear difference operator and the denominators of Ri, R 2 are d(r, s)g\{r, s) 
and d(r, s)g 2 (r, s), respectively. The recursion is the same as that obtained 
by applying Zeilberger's algorithm to the right hand side of ([5.5)1 . 
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